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]£2 Abstract 

Using pp collision data corresponding to 1.0 fb _1 integrated luminosity collected 
by the LHCb detector, two narrow states are observed in the yl^7r + 7r~ spectrum 
^ with masses 5911.97±0.12(stat)±0.02(syst)±0.66(ylg mass) MeV/c 2 and 5919. 77 ± 

0.08(stat)±0.02(syst)±0.66(yl£ mass) MeV/c 2 . The significances of the observations 
are 5.2 and 10.2 standard deviations, respectively. These states are interpreted as 
the orbitally excited A° b baryons, yl* (5912) and AJ°(5920). 
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i The system of baryons containing a b quark (beauty baryons) remains largely unex- 
plored, despite recent progress made at the experiments at the Tevatron. In addition to 
the ground state, A®, the baryon with the quark content bsd has been observed by 
the DO [l] and CDF [2] collaborations, followed by the observation of the doubly-strange 

5 baryon (bss) (3||4j^ The last ground state of beauty-strange content, (bsu), has 
been observed by CDF [5] . Recently, the CMS collaboration has found the corresponding 
? excited state, most likely with J p = 3/2+ |6|. Beauty baryons with two light quarks 

8 (bqq, where q = u,d), other than the A®, have been studied so far by CDF only Of the 

9 triplets with spin J = 1/2 and Z^' with J = 3/2 predicted by theory, only the 

10 charged states E^*^ have so far been observed via their decay to A^n^ 1 final states (7,8]. 

ii None of the quantum numbers of beauty baryons have been measured. 

12 The quark model predicts the existence of two orbitally excited A® states, Al°, with 

13 the quantum numbers J p = l/2~ and 3/2", respectively, that should decay to A®7t + it~ or 
u A^ry. These states have not previously been established experimentally. The properties of 

excited A® baryons are discussed in Refs. [9 15 . Most predictions give masses above the 



16 A®it + tt~ threshold, but below the E^it threshold. Observation of A* b Q states and measure- 

17 ment of their quantum numbers would provide a further confirmation of the validity of 
is the quark model, and the precise measurement of their masses would test the applicability 
is of various theoretical models used to describe the interaction of heavy quarks. 

20 This Letter reports the first observation of the Af states decaying into A®7i + 7i~ , and 

21 the measurement of their masses and upper limits on their natural widths. The data 

22 set of 1.0 fb _1 collected in pp collisions at the LHC collider at the center-of-mass energy 

23 = 7 TeV in 2011 is used for the analysis. 

24 The LHCb detector |16| is a single-arm forward spectrometer covering the pseudo- 

25 rapidity range 2 < r\ < 5, designed for the study of particles containing b or c quarks. 

26 The detector includes a high precision tracking system consisting of a silicon-strip vertex 

27 detector surrounding the pp interaction region, a large-area silicon-strip detector located 

28 upstream of a dipole magnet with a bending power of about 4Tm, and three stations 

29 of silicon-strip detectors and straw drift tubes placed downstream. The combined track- 

30 ing system has a momentum resolution Ap/p that varies from 0.4% at 5GeV/c to 0.6% 

31 at lOOGeV/c, and an impact parameter (IP) resolution of 20 \xm for tracks with high 

32 transverse momentum. Charged hadrons are identified using two ring-imaging Cherenkov 

33 (RICH) detectors. Photon, electron and hadron candidates are identified by a calorime- 

34 ter system consisting of scintillating-pad and preshower detectors, an electromagnetic 

35 calorimeter and a hadronic calorimeter. Muons are identified by a muon system com- 

36 posed of alternating layers of iron and multiwire proportional chambers. 

37 The online event selection (trigger) consists of a hardware stage, based on information 

38 from the calorimeter and muon systems, followed by a software stage which applies full 

39 event reconstruction. The software trigger used in this analysis requires a two-, three- or 
four-track secondary vertex with a high sum of the momenta transverse to the beam axis, 
Pt, of the tracks, and significant displacement from the primary interaction vertex (PV). 
In addition, the secondary vertex should have at least one track with p T > 1.7GeV/c, IP 
X 2 with respect to any PV greater than 16 (where the IP \ 2 is defined as the difference 



1 



44 of the PV fit y 2 with and without the track included), and a track fit xV n df < 2 where 

45 ndf is the number of degrees of freedom in the fit. A multivariate algorithm is used for 

46 the identification of the secondary vertices [TF| . 

47 The A® candidates are reconstructed in the A® — > A+tt~ , yl+ — > pK~7t + decay chain 

48 (addition of charge-conjugate states is implied throughout this Letter). The selection of 

49 A® candidates is performed in two stages. First, a loose preselection of events containing 

50 beauty hadron candidates decaying to charm hadron candidates is performed. It requires 

51 that the tracks forming the candidate, as well as the beauty and charm vertices, have 

52 good quality and are well separated from any PV, and the invariant masses of the beauty 

53 and charm candidates are consistent with the masses of the corresponding particles. 

54 The final selection requires that all the tracks forming the A® candidate have an IP 

55 y 2 with respect to any PV greater than 9, and the IP x 2 °f the A® candidate to the 

56 best PV (PV having the minimum IP x 2 f° r the A® candidate) is less than 16. Particle 

57 identification (PID) information from the RICH detectors is used to identify kaons and 

58 protons in the final state in the form of differences of logarithms of likelihoods between the 

59 proton and pion (DLL p7r ) and kaon and pion (DLL^) hypotheses. No PID requirements 
eo are applied to the pions from A® —> A+n~ decays to increase the A® yield: a significant 
ei fraction of these pions have momenta above 100 GeV/c where the PID performance is 

62 reduced. Finally, a kinematic fit is used which constrains the decay products of the A® 

63 and A£ baryons to originate from common vertices, the A® to originate from the PV and 



64 the invariant mass of the Aj candidate to be equal to the established A+ mass 18 



65 A momentum scale correction is applied to all invariant mass spectra in this analysis 



ee to improve the mass measurement using the procedure similar to 19 . The momentum 

e? scale has been calibrated using J/ip — > /i + \i~ decays, and its accuracy has been quantified 

es with other two-body resonance decays (T(1S) — > K® — > , — > K + K~). 

69 Signal and background distributions are studied using simulation. Proton-proton colli- 

70 sions are generated using Pythia 6.4 [20) with a specific LHCb configuration [21] . Decays 



7i of hadronic particles are described by EvtGen 22 in which final state radiation is gen 



72 erated using Photos 23 . The interaction of the generated particles with the detector 



73 and its response are implemented using the Geant4 toolkit 24 as described in Ref. 25 

74 The distribution of the A^ir~ invariant mass after the kinematic fit is shown in Fig. ]l| 

75 where a requirement of good quality of the kinematic fit is applied. In addition to the A® — >• 

76 A^tc~ signal contribution, the spectrum contains backgrounds from random combinations 

77 of tracks (random background), from partially-reconstructed decays where one or more 

78 particles are not reconstructed, and from A® — > A+K~ decays with the kaon reconstructed 

79 under the pion mass hypothesis. A fit of the spectrum yields 70 540 ± 330 signal events, 
so and the signal-to-background ratio in a ±25 MeV/ c 2 interval around the nominal A® mass 
si is S/B — 11. The fit to the A+n~ spectrum is only used to estimate the A® yield and the 

82 A® — > A+K~ contribution, and is not used in the subsequent analysis. 

83 The A® candidates obtained with the above selection are combined with two tracks un- 

84 der the pion mass hypothesis (referred to as slow pions from now on) to search for excited 

85 A® states. The tracks are required to have transverse momentum > 150MeV/c, and no 
se PID requirements are applied. A kinematic fit is applied that, in addition to all constraints 
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Figure 1: Invariant mass spectrum of A^tt~ combinations. The points with error bars are 
the data, and the fitted A® — > A+tt~ signal and three background components (A® — > A+K~ , 
partially-reconstructed and random background) are shown with different fill styles. 



described above for A® candidates, constrains the two slow pion tracks to originate from 
the PV and the invariant mass of the A® candidate to a fixed value of 5619.37 MeV/c 2 , 
which is a combination of the world average 18 and the LHCb measurement [26). The 
uncertainty on the combined A® mass obtained in this way, 0.69 MeV/c 2 , is treated as a 
systematic effect. Combinations with a good quality of kinematic fit, xV n df < 3.3, are 
retained. From the simulation study, this requirement is optimal for the observation of a 
narrow state near the kinematic threshold with signal-to-background ratio around one. 

The fit of the A^ir~ mass spectrum (Fig. [T| indicates the presence of the background 
from A® —> A+K~ decays at a rate around 12%, relative to the A® — > A+n~ signal. 
Alternatively, its rate can be estimated from the ratio of B + — > D°K + and B + — > D°ir + 



decays that equals to 8% 18 . Due to the A® mass constraint in the kinematic fit, the 
A®7r + 7i~ invariant mass distribution for this mode is biased by less than 0.1 MeV/c 2 if 
reconstructed under the A+n~ mass hypothesis, and has a resolution only a factor of two 
worse than that with the A^tt~ signal. After the kinematic fit quality requirement, the 
fraction of /i°7r + 7r~ with A® — > A+K~ decays compared to those with the A+tt~ is reduced 
to 8%. This mode is thus not treated separately, and its effect is taken into account as a 
part of the systematic uncertainty due to the signal shape. 

Combinations of A® candidates with both opposite-sign and same-sign slow pions are 
selected in data. The latter are used to constrain the background shape coming from 
random combinations of A® baryon and two tracks. The assumption that the shape of the 
background in A®7r + 7r~ and yl°7r ± 7r ± modes is the same is validated with simulation. The 
A®ii + tc~ and yl°7r ± 7r ± invariant mass spectra are shown in Fig. two narrow structures 
with masses around 5912 and 5920 MeV/c 2 are evident in the ylg7r + 7r~ spectrum. They 
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no are interpreted as the orbitally excited A® states, and are denoted hereafter as yl£°(5912) 

m and Af (5920). 

112 A combined unbinned fit of the A®ty + 'k~ and A^n^Ti^ 1 samples is performed to ex- 

n3 tract the masses and event yields of the two states. The background is described with a 

lu quadratic polynomial function with common parameters for both samples except for an 

us overall normalization. The probability density function (PDF) for each of the yl^°(5912) 

lie and yl^ (5920) signals is a sum of two Gaussian PDFs with the same mean. The rela- 

117 tive normalization of the two Gaussian PDFs are fixed to the values obtained from the 

us simulation of states with masses 5912 and 5920 MeV/c 2 and zero natural widths, while 

n9 the mean value and overall normalization for each signal are left free in the fit. The core 

120 resolution (width of the narrower Gaussian PDF) obtained from simulation is 0.19 and 

121 0.27 MeV/c 2 for A? (5912) and Af (5920), respectively. Study of several high-statistics 

122 samples (A®— > A^7r~, '0(25') — > J/V>/t + ?t~ , D* + — > D tt + ) shows that the invariant mass 

123 resolution in data is typically worse by 20% than in the simulation. Thus the nominal data 

124 fit uses the widths of Gaussian PDFs from the simulation multiplied by 1.2. The data 

12 5 fit yields 17.6 ± 4.8 events with mass M 4 *o (5912) = 5911.97 ± 0.12 MeV/c 2 and 52.5 ± 8.1 

126 events with mass M a *o/ 592 q) = 5919.77 ±0.08 MeV/c 2 . 



147 



127 Limits on natural widths T of the two states are obtained by performing an alternative 

128 fit where the signal PDFs are convolved with relativistic Breit-Wigner distributions. The 

129 dependence of Breit-Wigner width T on the A°/T + 7r~ invariant mass M is taken into 

130 account as T A *o(M) = T A *o x (q/qo) 2 x (M A *o/M). Here M A *o is the mass of the A* h ° 

131 state, and q^ is the kinematic energy for the decay of the state with mass M( A *oy. q^ = 

132 M(yj*o-) — M A o — 2M 7T , where M A o and M n are the masses of A® and 7r + , respectively. 

133 Scans of Breit-Wigner widths r yl .o( 5912 ) and T A *o^ 5920 ) are performed with all the other 

134 parameters free to vary in the fit. The upper limits are obtained without applying the 

135 mass resolution scaling factor of 1.2 as in the nominal fit to account for the uncertainty 

136 of this quantity: this gives a more conservative value for the upper limit. The 90% 

137 (95%) confidence level (CL) upper limit on T, which corresponds to 1.28 (1.64) standard 

138 deviations, is obtained as the value of T where the negative logarithm of the likelihood 

139 is 1.28 2 /2 = 0.82 (1.64 2 /2 = 1.34) greater than at its minimum. The 90% (95%) CL 

140 upper limit is r^»o (5912) < 0.66 MeV (0.83 MeV) for the ^°(5912) state, and r 4 *o (5920) < 

141 0.63 MeV (0.75 MeV) for the Af (5920) state. 

142 The invariant mass of the two pions, M(/T + 7r~), in the yl^°(5920) — > A®tt + tt~ decay 



143 is shown in Fig. |3j The background is subtracted using the s Weights procedure 27 

144 The weights are calculated from the fit to A b n + 7r~ invariant mass distribution, which is 

145 practically uncorrelated with M(vt + 7t _ ). The M(?t + /t _ ) distribution is consistent with 

146 the result of phase-space decay simulation, with x 2 / n df = 1.6 for ndf = 9. No peaking 
structures are evident. 



148 Systematic uncertainties on the mass measurement are shown in Table [T] The 

149 dominant uncertainty in the absolute A* b ° mass measurement comes from the uncer- 

150 tainty on the A® mass 5M A o = 0.69 MeV/c 2 ; it is propagated to the A* b ° mass un- 

151 certainty as 5M A *o = 5M A o x (M A o/M A *o) ~ 0.66 MeV/c 2 . This uncertainty mostly 
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Figure 2: Invariant mass spectrum of (a) A®tt + it~ and (b) A^it^tt^ combinations. The points 
with error bars are the data, the solid line is the fit result, the dashed line is the background 
contribution. 



cancels in the mass difference AM A *o = M A *o — M A o, where the residual uncertainty is 
5AM A *o = 5M A o x (AM A o/M A *o). The uncertainty of the signal parameterization is es- 
timated by using the simulated signal parametrization without applying the resolution 
scaling factor, by using the natural width for both states when left free in the fit, and by 
conservatively including the A® — > A+K~ contribution with the rate 12% parameterized 
from simulation. The uncertainty due to the background parameterization is estimated 
by: 

• using an alternative fit model for background description, 

• using the fit without the yl°7r ± 7r ± constraint, 

• using the fit with the background obtained from the simulation, 
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Figure 3: Invariant mass of the two pions from yl^°(5920) — > A®ir + ir~ decay. The points with 
the error bars are background-subtracted data, the solid histogram is the result of phase-space 
decay simulation. 



Table 1: Systematic uncertainties on the mass difference AM A *o between Al° and A®. 



Source of Systematic bias, MeV/c 2 

uncertainty AM A *o {5912) AM A *o (5920) 



A® mass 


0.034 


0.035 


Signal PDF 


0.021 


0.011 


Background PDF 


0.002 


0.002 


Momentum scale 


0.008 


0.013 


Total 


0.041 


0.039 



162 • fitting in the reduced invariant mass range 5910-5930 MeV/c 2 , 

163 and taking the largest difference from the nominal fit result as systematic uncertainty. 

164 The effect of the momentum scale correction is evaluated by varying the scale coefficient 

165 by its relative uncertainty 5 x 10~ 4 in simulated signal samples. 

lee The significance of the observation of the two states is evaluated with simulated 

167 pseudo-experiments. A large number of background-only invariant mass distributions 

168 are simulated with parameters equal to the fit result, and each distribution is fitted with 
leg models that include background only, as well as background and signal. The mean mass 

170 value of the signal PDF is not constrained in the fit to account for a trial factor in the 

171 range 5900-5950 MeV/c 2 . The significance is calculated as the fraction of samples where 
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172 the difference of the logarithms of fit likelihoods A log C with and without the signal 

173 is larger than in data. The fraction is obtained by an exponential extrapolation of the 

174 A log C distribution [28] that allows a limited number of pseudo-experiments to be used 

175 for a signal with high significance. The significance is then expressed in terms of the 

176 number of standard deviations (er). The significance of the yl^ (5912) state obtained in 

177 this way is 5.4<r for the A log C obtained from the nominal fit. To account for systematic 

178 effects, the minimum Alog£ among all systematic variations is taken; in that case the 

179 significance reduces to 5.2a. Similarly, the statistical significance of the ylj°(5920) state 

180 is 11.7(7, and the significance including systematic uncertainties is 10.2cx. 

181 The fit biases and the validity of the statistical uncertainties are checked with pseudo- 

182 experiments where the PDF contains both signal and background components. The fit 

183 does not introduce any noticeable bias on the measurement of the masses. The mass 

184 uncertainty for yi^°(5920) state is estimated correctly within 1% precision; however, the 

185 mass uncertainty for the yl^°(5912) is underestimated by 4%. This factor is taken into 
lee account in the final result. 

187 In summary, we report the observation of two narrow states in the A®7t + tt~ mass 

188 spectrum, yl^°(5912) and yl^°(5920), with masses 



200 
201 
202 



204 
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^4*°(59i2) = 5911.97 ± 0.12 ± 0.02 ± 0.66 MeV/c 
M A; o {5920) = 5919.77 ± 0.08 ± 0.02 ± 0.66 MeV/c 2 , 

189 where the first uncertainty is statistical, the second is systematic, and the third is the 

190 uncertainty due to knowledge of the A® mass. The values of the mass differences with 

191 respect to the A® mass, where most of the last uncertainty cancels, and the remaining 

192 part is included in the systematic uncertainty, are 

AilCi*o (5912) = 292.60 ± 0.12(stat) ± 0.04(syst) MeV/c 2 , 
AM^op) = 300.40 ± 0.08(stat) ± 0.04(syst) MeV/c 2 . 

193 The signal yield for the /lj* (5912) state is 17.6 ± 4.8 events, and the significance of 

194 the signal (including systematic uncertainty and trial factor in the mass range 5900- 

195 5950 MeV/c 2 ) is 5.2 standard deviations. For the Af(5920) state, the yield is 52.5 ± 8.1 

196 events and the significance is 10.2 standard deviations. The limits on the natural widths of 

197 these states are r yl ,o (5912) < 0.66 MeV (< 0.83 MeV) and r 4 *o (5920) < 0.63 MeV (< 0.75) 

198 at the 90% (95%) CL. 

199 The masses of Al° states obtained in our analysis are 30-40 MeV/c 2 higher than in 



the prediction using the constituent quark model 12 , and 20-30 MeV/c lower than the 
predictions based on the relativistic quark model |11| , modeling the color hyperfine inter- 



action 14 and an approach based on the heavy quark effective theory 15 . Calculation 
203 involving a combined heavy quark and large number of colors expansion (9||10) gives a 
value roughly in agreement, although only the spin-averaged prediction is available. The 



205 earlier prediction based on the relativized quark potential model 13 matches well the 

206 absolute mass values for both states, but the A® mass prediction using this model is 

207 35 MeV/c 2 lower than the measured value. 
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